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ABSTRACT 

The mass of molecular gas in an interstellar cloud is often measured using line 
emission from low rotational levels of CO, which are sensitive to the CO mass, and 
then scaling to the assumed molecular hydrogen H2 mass. However, a significant 
H2 mass may lie outside the CO region, in the outer regions of the molecular 
cloud where the gas phase carbon resides in C or C^. Here, H2 self-shields or 
is shielded by dust from UV photodissociation, whereas CO is photodissociated. 
This H2 gas is "dark" in molecular transitions because of the absence of CO and 
other trace molecules, and because H2 emits so weakly at temperatures 10 K 
< T < 100 K typical of this molecular component. This component has been 
indirectly observed through other tracers of mass such as gamma rays produced in 
cosmic ray collisions with the gas and far-infrared/sub millimeter wavelength dust 
continuum radiation. In this paper we theoretically model this dark mass and 
find that the fraction of the molecular mass in this dark component is remarkably 
constant (~ 0.3 for average visual extinction through the cloud Ay — 8) and 
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insensitive to the incident ultraviolet radiation field strength, the internal density 
distribution, and the mass of the molecular cloud as long as Ay, or equivalently, 
the product of the average hydrogen nucleus column and the metallicity through 
the cloud, is constant. We also find that the dark mass fraction increases with 
decreasing Ay, since relatively more molecular H2 material lies outside the CO 
region in this case. 

Subject headings: ISM: clouds 



1. INTRODUCTION 



Various observations have indicated that a substantial amount of interstellar gas exists in 
the form of molecular hydrogen (H2) along with ionized carbon (C+), but little or no carbon 
monoxide (CO). The total mass in molecular hy drogen has been estimated from gamma ray 



observations from C OS-B (iBloemen et al.l 119861) and EGRET (Energetic Gamma Ray Ex- 
periment Telescope) ( Strong fc Matto:?3 19961 ) and analysis of this data showed more gas mass 



than can be accounted for in H I and CO alone ( IGrenier et al.ll2005[ ). In addition, the dust 



column density maps of the Galaxy from DIRBE, and maps of the 2MASS J-K extinction 



show additional gas not seen in H I or CO ( iGrenier et al.ll2005[ ) and is presumably molecular 



hydrogen. In molecular line observations and modeling of low column density molecular 
clouds the H2/CO ratio is found to be variable and much larger than 1 0'^ so that only a 
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clouds along with H I and CO observations, found H2 masses comparable to the H I masses 
with only small amounts of CO. A mixture of H2 and C+ is also inferred to exist in dif- 
fuse clouds where the H2 and CO columns are me asured by UV absorption spectroscopy 



and the CO accounts for only a trace amount of C (ISonnentrucker et al.l 120071 : iBurgh et al. 



20071 : ISheffer et al.l 120081). T his component of the interstellar medium (ISM), termed "dark 
gas" by iGrenier et al.l (120051 ) . is also inferred to exist in extragalactic observati ons compar- 
ing far-infrared a nd CO mass estimates especially in low metallicity galaxies ( llsraell Il997l : 
Lerov et al.ll2007l ). 



Such an H2 and C^ la yer is also predicted from theor etical models of diffuse gas 



( Ivan Dishoeck fc Blacklll988l ) and surfaces of molecular clouds (ITielens fc Hollenbachlll985l ) 
that indicate the transition from C^ to CO is deeper into the cloud than the transition 
from H to H2. Essentially, the theoretical models show that H2 self-shields itself from UV 
photodissociation more effectively than CO. This layer is "dark" in rotational H2 transitions 
primarily because the ground state transition 0-0 S(0) at 28 /xm lies about I\Ejk ~ 512 K 
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above ground; at the temperatures 10 K < T < 100 K typical of the dark component the 
fluxes from the H2 rotational transitions lie below current sensitivities. This layer is dark in 
CO because of its very low abundance. Although the dark gas layer is "dark" in H2 and CO, 
it does emit mainly in [C II] 158 /xm fine-structure line emission. In preliminary estimates 
(M. Wolfire et al. 2010, in preparation) the calculated local Galactic [C II] emission from the 
WIM and CNM (diffuse phases that are not associated with molecular clouds) underestimate 
the COBE observations of the line emission in the plane by a fact or of 1/3 to 1/2 and could 



be another indicator of dark gas. Similarly, IShibai et al.l ( 1l99ll ) and ICubick et al.l (120081 ) 



found the bulk of the [C II] emission in the Galactic plane seen by BIRT and COBE arises 
in neutral gas associated with molecular clouds. 

In this paper we present models of molecular cloud surfaces to estimate the mass of gas 
in the "dark" component. In §2 we discuss modifications to existing photodissociation region 
(PDR) codes, and the modeling procedure. In §3 we define the dark gas mass fraction and 
other parameters used in our models, discuss the average gas density distribution assumed 
in our modeling, and derive an expression for the dark gas fraction in terms of parameters 
found in our numerical modeling procedure. The modeling results start in §4 with a simple 
isobaric cloud model in the limit of clumps that are optically thin to the incident radiation 
field, resulting in an estimate of the dark gas fraction for a typical giant molecular cloud. 
We then in §5 enhance this model with the inclusion of turbulent pressure and find the dark 
gas fraction as a function of incident field strength and cloud mass. We also discuss the 
variation in the dark gas fraction as a function of metallicity (over a limited range), the 
average hydrogen column through a cloud, the average visual extinction through a cloud, 
and the opacity of the clumps. We compare our results with observations. Finally, in §6 we 
conclude with a discussion and summary. 



2. Models 

Our goal is to determine the mass in the C+/H2 layer — the dark gas — in molecular 
clouds. Our analysis is based in part on the r e sults of a modified version of the photodisso- 



ciation region (PDR) code of iKaufman et al.l ( l2006l ). In applying this code, we are making 



four main approximations: First, we assume that the inhomogeneities in actual molecular 
clouds can be approximated as clumps of density ric occupying a fraction /v of the volume; 
we neglect the interclump medium, so that the locally volume-averaged density, n is given 
hy n = fyUc. Second, we assume that these clumps are optically thin to UV radiation; 
as a result, the radiative transfer is the same as that in a homogeneous medium with the 
mean density, n. The validity of this assumption is discussed in §5.41 below, where we show 
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that even in the limit of very optically thick clumps, our models give essentially the same 
dark gas mass fraction. Third, in the PDR solution to the variation with depth of chemical 
abundances and gas temperature, we assume that the layer of dark gas is not geometrically 
thick, so that we can apply the results of a ID slab model to the spherical problem. However, 
in calculating the mass of gas in the dark gas layer, we do take into account the spherical 
geometry. Fourth, we assume that we can calculate the chemical abundances in steady state. 
This assumption is discussed in §5.51 after we have derived densities, turbulent speeds and 
characteristic distances. With these approximations, the physical conditions and chemical 
abundances within the cloud are a function of the optical depth from the cloud surface and 
we can apply a single, one-dimensional PDR model to a continuous density distribution in 
which the gas density ric inside the clumps can vary as a function of the distance r of a 
clump from the center of the GMC 

The PDR models find the gas temperature in thermal balance and the chemical equi- 
librium abundance of all dominant atomic and molecular species as a function of depth into 
a gas layer. The models require the incident far-ultraviolet (FUV; 6 eV < hv < 13.6 eV) 
radiation, extreme ultraviolet (EUV, 13.6 eV < hv < 100 eV) radiation, soft X-ray (100 eV 
< /iz/ < 1 keV) radiation, cosmic-ray flux, and either the spatial distribution in density or 
the spatial distribution in thermal pressure. If the density distribution is provided, the code 
will iterate on the gas temperature at fixed density until thermal balance is achieved. If the 
thermal pressure {Pth/k = rifT, where rit is the total particle density) is provided, the code 
will iterate on both the density and temperature. With these inputs, the temperature, and 
the atomic and molecular abundances are calculated self-consistently. 



2.1. Modified PDR Code 



The PD R code of iKaufman et al.l (20061) is based on two main parts. The first part is 
derived from Kaufman et al. ( 19991 ). Wolfire et al. (1990), and Tielens fc HoUenbach ( 1985 ) 
models and is the main program for calc ulating the chenaistry, thermal balance, and line 
emission. The second part is based on the iLe Petit et al.l (120061) Meudon PD R code and is 
used for all of the molecular hydrogen processes. The iKaufman et al.l (120061 ) code assumes 
a ID UV flux normal to the cloud surface. Here we consider the average radiation field 
produced by all OB associations illuminating the giant molecular cloud (GMC) and the field 
at the cloud surface is expected to be relatively isotropic over 27r steradians. In order to 
account for the isotropic field, we assume a ID flux incident at an angle of 60 degrees to 
the normal. This unidirectional field has the same flux normal to the cloud surface as the 
isotropic field, and thus the same energy input to the cloud, and it more closely approximates 
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the depth penetration of the field. The photo rates are modified to account for twice the path 
length of the incident field for a given Ay measured normal to the surface. In the Meudon 
portion of the code, the isotropic field is accounted for explicitly in the depth dependence of 
the photo rates. 

We have also included several changes in the elemental abundances and chemical network 
that affect the ca r bon c hemistry. We adopt the gas phase carbon abundance = 1-6 x 10"*^ 
from ISofia et al.l ( 120041 ). where = nc/n and where nc is the gas phase number density 
of elemental carbon and n is the number density of H nuclei. W e also include the C"*" 
radiative and dielectronic recombination rates from Badnell ( 2006a ) and (Badnell 2OO6b)0, 
who have found that t he rates at low ternperat ures can be a factor of 2-3 higher than was 



used in the past (e.g., Nahar fc Pradhan 1997h. Another change in the chemical network 



is to include OH production on grains from iHoUenbach et al.l (|2009[ ). We adopt the fastest 
possible rate, having each collision of O with a grain combine with H and produce OH. This 
limit tends to maximize CO production and minimize the dark gas fraction. We also include 
th e ion-dipole reaction rat es for C^, , and He"*", reacting with OH and H2O as discussed 



m 



Hollenbach et al.l (120091). 



2.2. External Radiation Field 

The external fiuxes of FUV radiation, EUV/soft X-ray radiation, and cosmic rays are 
model inputs and must be specified. Associations produce bright photodissociation regions 
on nearby cloud surfaces. For example, the Trapezium cluster in Orion produces the bright 
Orion PDR on cloud surfaces only ~ 0.1 pc distant. As an association evolves, the Stromgren 
region expands, the embedded H II region breaks out of the cloud, the resultant champagne 
fiow and blister evaporation disperses the cloud and the distance between the cloud and 
association increases, thus allowing more of the cloud to be illuminated than the initial hot 
spot. The combination of cloud and association evolution, the distribution of associations 
around the cloud, and the continuous formation of new associations, elevates the average 
radiation field illuminating the cloud to a field strength that is greater than the average 
interstellar radiation field. A detailed evaluation of this process is presented in Wolfire, 
Hollenbach, & McKee (2010, in preparation). They find the average field incident on clouds 
is o n the order o f ~ 20 times higher than the local Galactic interstellar fi eld. We adopt 



the iDraind (119781 ) interstellar radiation field that is equal to ~ 1.7 times the iHabingl (119681 ) 



local interstellar field (1.6 x 10 ^ erg cm ^ s ^ for FUV photons). We also use the notation 



^http://anidpp.phys. strath. ac.uk/tamoc/DR 
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that Go is the radiation field in unit s of th e Habing field and Gq is the field relative to the 
local Galactic field found by iDraind ( 119781 ). so that Gq = Gn /1.7. To clarify the de f initio n 
of Gq, which is consistent with the original definition given in lTielens fc HoUenbachl ( 119851 ) . 
Gq is the ratio of the photorates at the surface {Ay = 0) of our cloud to the same rates in 
the Draine field in free space. Thus, Gq = {2711)/ {An Id) = 0.51 /Id, where / is our incident 
isotropic intensity and Id = 2.2 x 10~^ erg cm^^ s^^ sr^^ is the Draine intensity for the local 
ISMo Here we examine a range of FUV fields, 3 < Gq < 30, consistent with (and generously 
encompassing) the results of Wolfire et al. (2010). Throughout this paper w e simply scale 
the fiux by Gq and keep the spectral energy distribution fixed as given by the iDraind (Il978[ ) 
field. 



As discussed in lWolfire et al.l (|2003[ ). the soft X-ray/EUV fiux from the diffuse interstel- 
lar field can be produced by a combination of supernova remnants and stellar sources. The 
field shining on the molecular cloud can be enhanced above the interstellar flux due to the 
OB associations near the molecular cloud. We assume that we are outside any of the H II 
regions associated with the OB associations, and simply scale the soft X-ray/EUV fiux by 
the same factor that we do for the effective far- ultraviolet radiation field Cxr = Gn. A corol- 



l ary o f this assumption is that the fraction of gas that is ionized is negligible. I Wolfire et al. 



( 120031 ) adopted a standard soft X-ray/EUV shielding column of neutral gas of column density 



N = Ix 10 cm~ . The layer N < 10 cm~ mainly consists of WNM gas photoevaporated 
from the molecular and CNM clump surfaces and does not provide significant shielding to 
the FUV photodissociating radiation field. We are mainly concerned with the deeper layers 
and assume that the clouds are shielded by a column of this order. As we shall demonstrate 
in ^ the absorption of X-ray/EUV radiation at columns in the range 10^^ cm~^ < ~ 10^^ 
cm~^ leads to a drop in the thermal pressure by a factor of ~ 10 in the H I layer in the 
cloud. 

The cosmic-ray fiux is also associated with massive star evolution, but is not as likely 
to be simply proportional to Gq as for the EUV/soft X-ray fiux. Since the temperatures of 
cloud interiors do not exceed T = 10 — 20 K, the cosmic ray fiux in cloud interiors cannot 
be higher than a factor of ~ 2 — 3 times the diffuse ISM value or else the heating rate 
would push internal cloud temperatures higher than observed. In most of our calculations, 
we therefore assume that cosmic-ray rates are not enhanced by the nearby molecular cloud 
OB associations. 



The cosmic ray rate can also affect the formation of CO in the diffuse gas and outer 



^Note that an opaque cloud in the local ISRF therefore experiences Gq = 0.5, since the molecule at the 
surface is only illuminated from 27r sterradians. 
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molecular cl oud layers through the production of OH as discussed by Ivan Dishoeck fc Black 
( ll986l . Il988[ ). There are several observations, mainly of Hj]", that indicate hig her cosmic ray 
ioniz ation rates by a factor of ~ 10 along select lines of sight in the diffuse ISM ( llndriolo et al. 



20071 ): however the B.^ observatio ns do not indicate higher cosmic ray rates in the denser, 



molecular clouds considered here ( iMcCall et al.lll999l ). We adopt our standard rate for the 
majority of this paper, but consider the effects of higher cosmic ray rates in the outer portions 
of the cloud in §5.2.6[ 



3. Density Distribution and Dark-Gas Fraction 



To estimate the dark-gas fraction in a molecular cloud, we assume that it is spherical. 
Let Rco be the radius of the CO part of the cloud, defined by the condition that the optical 
depth from Rqo to the surface in the J = 1 — transition is rco = 1- Let > Rco 
be the radius of the part of the cloud in which hydrogen is molecular; we define this as the 
radius at which the mass density in H2 molecules equals that in H atoms. Finally, the total 
cloud radius, -Rtot, includes the H I shielding layer that absorbs the FUV radiation incident 
upon the cloud, enabling the gas to become molecular. The corresponding masses are the 
CO cloud mass, M{Rco), the total molecular mass, M{Rii^), and the total cloud mass, 
^(-Rtot), respectively, including the associated helium. The dark gas extends from Rqq to 
the atomic- molecular transition radius at Ru2- The dark-gas mass fraction is then 

_ MjRn,) - MjRco) 

Note that /dg is the fraction of molecular (H2) gas (and not the fraction of the total mass, 
atomic plus molecular) that is dark. It is also not the ratio of the dark gas to the gas in the 
CO zone. 

The gas in the cloud generally has large density fluctuations. In the H I, this is because 
the gas can exist i n two distinct phase s, cold (T ~ 10^ K) and dense, and warm (T ~ 10^ K) 



and tenuo us fe.g.. IWolfire et al.ll2003[) . whereas the molecular gas is generally supersonically 



turbulent f Larson 198ll ) and therefore has large density variations. Let n(r) be the density of 



H nuclei at a point in the clou d, and l e t n(r) be the locally volume-averaged density at that 
radius in the cloud. Following iLarsoru ( 119811 ). we assume that n(r) oc 1/r. This corresponds 
to a global column density that is independent of radius, 

M(r) 



N 



2rn{r) 



const, 



(2) 



where /xh 

use a column density N 



2.34 X 10 g is the mass per H nucle us. For numeric a l eval uation, we usually 
1.5 X 10^2 



cm from ISolomon et al.l (Il987l ) . The column N 
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corresponds to an average visual extinction Ay through the cloud given by 

NZ' N 

= TO in2i 2 = = 5.26iV22Z', (3) 

1.9 X 10^^ cm-^ No 

where Nq = 1.9 x 10'^^Z'~^ cm~^ is the radial column den sity corresponding to a visual 
extinction of unity jBohlin et al.lll978[ iRachford et"aD 120021 ) . N22 = N/iW^^ cm'^), Z' = 



Z/Zq is the metallicity relative to solar, and assuming that the gas to dust ratio is constant 
within the cloud and the dust opacity scales with the gas metallicity Z'. The relation between 
the mass and the locally volume-averaged density is 



^3/2 

n(r) = 30.4 cm'^ (4) 

where Meir) = M(r)/(10^ Mq). 

The PDR code is run for a cloud of mass M(i?co) on a grid of visual extinction, Ay, 
measured from the surface of the cloud at i?tot to a point r inside the cloud. We denote the 
corresponding column density as 



N{r) = - ndr' (5) 
J Rtot 

to distinguish it from the globally averaged value in equation (|2]). Note that N{r) is the 
column density outside of r, whereas N is the average column density through the cloud 
and does not depend on r for our 1/r average density distribution. The visual extinction 
Ay corresponding to a column density N is given as in equation (|3]) with N replaced by A^. 
Note that the column density, (or extinction. Ay) between i?co and -Rtot is not known a 
priori and is one of the main results of our calculation. In addition, for a 1/r mean density 
distribution, the column through the cloud center is infinite, but the total mass within the 
Tco = 1 surface is well defined and is one of the input parameters. In practice we run the 
PDR model to a depth of Ay = 10 and adjust the inner radius of our grid at Ay = 10 so that 
the integrated CO cloud mass, M(i?co), is the input value. Our results do not depend on 
the value of Ay we choose for the inner radius since there is so little mass or radius contained 
in the small central region where the radial Ay goes to infinity as r goes to zero. 

The results from the PDR code determine both the extinction at the H I - H2 interface, 
Ay{Rn^) and the extinction at the surface of the CO cloud, Ay{Rco)- The difference between 
these gives the extinction of the dark-gas layer. 



(6) 
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For a cloud with an r ^ density profile, corresponding to M(r) oc r^, the dark-gas mass 
fraction is 

/DC = (7) 



R 



Ha 



1 - exp ( ) , (8) 

^7T> (9) 



Z'N, 



22 



= l-exp(^ ^- J. (10) 

This result is generalized to other density profiles in Appendix [Aj The dark-gas fraction 
thus depends upon only two quantities, the radial column density of the dark-gas layer, 
AiV = iVoAAy,DG; and the mean column density of the cloud, N . When expressed in terms 
of the extinction of the dark gas layer, AAy dg; a dependence on the metallicity, Z' , also 
enters. Perhaps most simply, when expressed in terms of AAy ^G and Ay-, the fraction 
depends only on the ratio AAv^dg/^v- This is intuitive since (discussed later in detail), 
AAy^DG is a measure of the dark gas mass and Ay is a measure of the molecular mass. 



4. Radiatively Heated Clumps at Constant Pressure 

As a first step, we consider the case of isobaric clumps with a thermal pressure that 
is independent of position within the cloud. The clumps are heated by an FUV/soft X-ray 
radiation field that is attenuated by the gas between the clumps and the surface of the 
PDR, plus a column of 10^^ cm~^. (The clumps are also heated by the cosmic rays, but 
this heating is generally less important in the dark-gas region.) The thermal pressure is 
-fth = XtUckT, where Xt is the sum over the abundances of all species relative to hydrogen 
nuclei, Xt = TiiUi/n] for atomic g ~ 1.1, while for molecular g ^ 0.6. Note that 
the density derived from the isobaric PDR model is ric, the density of the cold (T ^ 300 
K) gas component. We assume that warm (T ~ 8000 K) H I at the same thermal pressure 
fills the remaining space, but we assume that this has a negligible fraction of the mass. As 
noted in ^ above, the volume filling factor is given by /v(r) = n{r)/nc{r), where n(r) is the 
locally volume-averaged H nucleus density (i.e., the density at r averaged over the clumps 
since we are neglecting the interclump medium) and nc{r) is the H- nucleus density of the 
cold (clump) gas. 

The model outputs include the local density in the clumps, nc{r), the volume filling 
factor of the clumps, /v, the fraction of atomic hydrogen, xhi = nm/nc, and the fraction of 
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molecular hydrogen, xn^ = nn^/ric = |(1 — xhi) as functions of extinction from the surface 
of the cloud, Ay- The total mass in each component is found by integrating the density 
distributions 

/•Rtot 

M{Rn,)= 2/iHXH,(r')ne(r')/v(r')47rr'2dr' , (11) 

^0 

pRco 

M{Rco)= 2fiuXu,{r')n,{r')fy{r')A7Tr'^dr' , (12) 

^0 

where Rqo is the radius of the tco = 1 surface, and Rtot is the outer radius of the entire 
cloud, which includes (from center outwards) the CO region, the region with H2 and C"*", 
and the outer atomic envelope with atomic H and C+ (see Fig. [T] for an illustration of the 
various radii and optical depths). Note that -Rco is an input to our model while i?tot is an 
output; it is the outer radius that gives just enough shielding such that CO J=l-0 becomes 
optically thick at -Rco- 

We have run cases for a representative cloud of mass M{Rqo) = 1 x 10® M©, with 
incident radiation fields G'q = Cxr = 10, and two fixed values of the pressure, Pth/k = 10^ K 
cm~^ and 10^ K cm~^, which covers the observed range of thermal p ressure deduced from 
^^CO and ^^CO observations of molecular clouds in the Galactic plane (ISanders et al.lll993l ). 



For Pth/k = 10^ K cm~^, the atomic-molecular transition occurs at ^^(i^Ha) = 0.54, the 
transition to optically thick CO (at rco = 1) occurs at Av{Rco) = 1-2, and the dark-gas 
fraction is found to be /do = 0.28. The total molecular mass is related to M(i?co) by 

M(fl„,) = , (13) 

-L - /DC 



with M(i?H2) = 1.4M(i?co) in this case. For Pth/A; = 10^ K cm'^, we find AviRn^) = 0-10, 
Av{Rco) = 0.86, /dm = 0.31 and M{Rn^) = 1.4M(i?co)- Although the transitions to H2 
and CO gas are shifted to lower column densities for the Pth/k = 10^ K cm~'^ case, the 
dark-gas faction remains relatively unchanged. 

To investigate the sensitivity of this result to variations in the density distribution and 
FUV radiation fields, we have constructed more realistic cloud models that include the effects 
of a two-phase thermal pressure and the effects of turbulence. 



5. Clumps in a Turbulent, Two-Phase Medium 

The illumination of the cloud by X-ray/EUV and FUV radiation will heat and evapo- 
rate the clumps, leading to a two-phase region in the outer regions of the molecular cloud 
consisting of warm T ~ 8000 K interclump gas and cooler T < 300 K clumps. We assume 
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the two-phase thermal pressure is the geometric mean of the minimum and maximum ther- 
mal pressures allowed for a two-phase medium, P^^ = (-Pmax-Pmin)^''^) and we again assume 
that the mass of the interclump gas is negligible. Shocks or turbulence will tend to keep 
the thermal pressure above Pmim while condensation of WNM into CNM will tend to keep 
the thermal pressure below Pmax- The geometric mean is roughly consistent with the mean 
thermal pre ssure in cold gas found i n numerical simulations of two-phase ISM disks with 
turbulence (jPiontek fc Ostrikerl 120071 ) . For sufficiently strong external radiation fields, pho- 
toevaporation of the cool clumps near the outer boundary of the cloud increases the density 
(and therefore the thermal pressure) of the all-pervasive interclump medium, so that the 
thermal pressure due to radiative heating exceeds the turbulent pressure there. However, 
the thermal pressure due to radiative heating drops from the cloud surface inward due to 
absorption of the soft X-ray/EUV and FUV radiation. 

For sufficiently weak radiation fields such as exist in the interior regions of molecular 
clouds, most of the gas is cold and the turbulent motions a re very superson i c, wh i ch generates 
a wid e range of therm al pressures. Both observational (jLombardi et al.l l2006l : iRidge et al. 
20061 ) and theoretical (I Vazquez- Semadenil Il994j : lOstriker et al.l 1200 ll ) studies suggest that 
turbulent clouds have a log -normal density distributio n that applies to the dense component 
in the two-phase medium (lAudit fc Hennebelldl2010l ). The mass- weighted median density, 
(^)med, in this turbulent medium is related to the volume-averaged density, n, by 



(?^)med = nexp(/i) 



(14) 



where 



/i^0.51n(l + 0.25M^ 



(15) 



( Padoan et al.lll997l ). Ai = y/Sa/cg is the 3D Mach number, a is the ID velocity dispersion, 
and Cs is the sound speed. The typical thermal pressure in a turbulent medium is P^^'^^ = 
Xt{n)j^e^kT: half the mass is at a greater thermal pressure than this and half at a lower one. 
To treat the complex situation in which both radiative heating and turbulence determine 
the pressure, we assume that the typical thermal pressure in the gas is the larger of that due 
to radiative heating, P^^, and that due to turbulence, P^^'^^, 



Pth = max (P,^ Pr') • 



(16) 



To determine the Mach number, we note that the relation between the linewidth and 



size of a cloud is given by the identity 



a{r) = 0.52 



102 Mq pc- 



r 

pc 



1/2 



km s 



(17) 
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where Ovir = ^o'^rjGM is the virial parame ter and S = Mli^r^ = i ii^N is the average 
mass surface density of the cloud (see eq. 27 in lMcKee &: Ostrikerll2007l ). iHeyer et ahl (120091 ) 
emphasized the importance of the dependence of the hnewidth on column density and showed 
that the data on Galactic mol ecular clouds a r e con sistent with this scaling with a^ir = 1- 
By themselves, the GMCs in ISolomon et al.l (Il987l )'s sample do not show clear evidence 
for the iV^/^ scaling; that b ecomes evident onl y when higher resolution data are included. 
Nonetheless, no rmalizing the iHeyer et al.l (120091 ) relation so that it agrees with the linewidth- 
size relation of Solomon et al.l ( 19871 ) at the mean column density of the latter 's sample 
[N = 1.5 X 10^^ cm~^) gives a linewidth-size relation 



air 



0.72 



22 



1.5 



1/2 



1/2 

— I km s~^ 



We adopt this relation for the ID velocity dispersion here. The corresponding 3D Mach 
number is = y/3a{r)/cs, where the isothermal sound speed is Cs(r) = [xtkT / {jju)]^^"^ ■ 

At each Ay step we use a precalculated lookup table to find the appropriate two- 
phase thermal pressure P^^ = (-Pmax-Pmin)^^^, as a function of total column density from 
the cloud surface and molecular fraction f(B.2) = 2N }^^/N (see Fig. |2]). The pressure is 



calculated in a manner similar to IWolfire et al.l ( 120031 ) except that here we include the H2 
self-shielding to depth Ay- The maximum, Pmax, and minimum, Pmm, thermal pressures for 
a two-phase medium are found from phase diagrams of thermal pressure versus density. In 
the regime where two-phase pressure dominates, the gas temperature and local density nc{r) 
are found self-consistently by iteration. In addition, however, we need the mass-weighted 
median density in the turbulent medium, (n)med (eq. [H]), which is a function of the cloud 
temperature T(r) [through Cs(T)], and the cloud radius [through 0"(r), n(r), and T{r)]. The 
solution requires two sets of iterations. First, we iterate between (n)niod and T at each Ay 
grid point to find a self-consistent density. Second, we convert the Ay grid to a radius grid 
and rerun the temperature, density iteration. 

We note that the H2 formation rate per unit volume is proportional to rimn, so one 
might expect that the mass- weighted mean density (n)M = (^^)v/?^ = nexp{2fi) would be 
appropriate for H2 formation and thermal balance rather than the median density, (n)med = 
nexp(/i); here the subscript "V" refers to the volume average. However, in highly supersonic 
gas, only a very small fraction of the mass is above the mass-weighted mean density, so the 
bulk of the chemistry occurs in the lower density gas. When the transition from atomic to 
molecular gas is sharp, then the conditions for the transition are determined by the lower 
de nsity atomic gas, ii o t by the molecular gas. We also note that numerical simulations 
of Glover fc Mac Low J20071 ) suggest that II2 formation proceeds rapidly in a turbulent 
medium where the H2 forms in high density gas and remains molecular when transported to 
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a lower density region. However, using {n)M for the typical density at which the formation 
occurs does not take into account the fact that in very dense regions the gas becomes fully 
molecular, so that the formation rate drops; thus use of the higher density would overestimate 
the formation rate. Therefore, we use {n)med for the typical density of a clump. 



5.1. Results for Clumps in a Turbulent Medium 

We have applied our optically thin clump models to molecular cloud masses M(i?co) = 
1 X 10^ Mq, 3 X 10^ Mq, 1 X 10^ Mq, and 3 x 10^ M© for incident radiation fields G'q = Cxr = 3, 
10, and 30. The lower end of the mass range is chosen to be the lowest cloud mass that 
could still produce OB stars sufficient to illuminate, heat, and dissociate the molecular cloud 
with FUV radiation . The upper end is about 1/2 the maximum mass cloud in the Galaxy 



( IWilliams fc McKed 119971 ). The range in radiation fields generously covers the expected 



elevated field from the sum of nearby stellar associations. 

We show in Figure [3] (top) the calculated thermal pressures for a representative cloud of 
mass M{Rco) = 1 x 10^ Mq with incident radiation fields G'q = Cxr = 10- this case, the 
two-phase thermal pressure (i.e., the thermal pressure due to radiative heating) dominates 
in the outer portion of the cloud, but drops initially due to absorption of EUV/soft X-rays 
and then due to absorption of FUV by d ust. We note t hat t he drop in two-phase pressure 



is greater than expected from Table 3 in IWolfire et al.l ( 120031 ). They quote a decrease by a 
factor ~ 1.8 for P^^ (given as Pave in that paper) between N = 10^^ cm~^ and = 10^*^ 
cm~^, whereas we find a drop by a factor of ~ 3.7. The difference is due to the treatment of 
cosmic rays: here we do not scale the cosmic ray ionization rate with G'q but hold it fixed at 
the local Galactic value . Thus the ionization falls off faster with depth into the cloud than in 



the lWolfire et al.l (|2003[ ) models. The turbulent thermal pressure dominates for Ay > 0.001. 
The transition to H2 (at n-a^/n = 0.25) occurs at Av{R]i2) = 0.34, the transition to optically 
thick CO (at tco = 1) occurs at Av{Rco) = 1-0, and the dark gas fraction is /dg = 0.28. 

The calculated densities are shown in Figure^ {bottom). We denote the clump density 
that would occur in two-phase equilibrium in the absence of turbulence by n2p. For the 
actual clump density, we take 



max(n2p, (n)med), (19) 



where the turbulent mass-weighted median density, (n)med, is given by equation f|T^ . For 
the case in the figure, we have Uc = for Ay < 0.001 and ric = {n)^ed for Ay ^ 0.001. 
Both H2 and CO form well within the turbulence-dominated region, a result that holds for 
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all cloud masses and radiation fields considered herejfl Note that the volume filling factor 
/v = n/uc ~ 0.1 in the outer part of the cloud where the thermal pressure dominates and 
two phases exist. Here, WNM fills the rest of the volume with a density of ~ 10~^?t,c. This 
implies the mass of the interclump medium is negligible (~ 0.1) compared to the mass in 
clumps, as we have assumed. In the turbulence-dominated region, there is no such thing 
as an interclump medium, but rather a distribution of densities that fill the volume. Here, 
most of the mass is at densities near (n)ined, as assumed. 

Figure H] also shows the temperature distribution {top) and the distribution of abun- 
dances for C"*", C°, CO, and OH [bottom). The CO amounts to only ~ 30% of the total 
carbon abundance when the J=l-0 line becomes optically thick at a CO column density 
of iV(CO) ^ 2 X 10^^ cm~^. Note tha t we do not include fr eeze out of H2O on grains, 
which would be important at Ay ^ 3 (jHoUenbach et al~ 2009 ). so our OH chemistry be- 



comes somewhat unreliable at these depths. However, this happens after CO is formed and 
will not change our results. 

For comparison with the M(i?co) = 1 x 10^ Mq case we also show in Figures and E] the 
results for M(i?co) = 1 x 10^ Mq. Figure E] shows the thermal pressures and densities and 
Figure [6] shows the gas temperature and abundances. The lower cloud mass results in higher 
n (eq. H]) and thus higher clump density (eq. [T^ and higher turbulent thermal pressure. The 
turbulent thermal pressure dominates the two-phase thermal pressure at all Ay. We find 
AviRn^) = 0.23, Av{Rco) = 0.95, and the dark gas fraction /dg = 0.30. 

The optical depths Ay(i?H2) and Av{Rco) are presented in Figure [7] (top) for cloud 
masses M{Rco) = 1 x 10^ M© and M(i?co) = 3 x 10^ Mq and for G'^ = Cxr = 3, 10, and 
30. The optical depths are measured from the cloud surface inward. Figure [7] ( bottom) shows 
the cloud radii (measured from the cloud center) corresponding to Ay(i?H2); ^v{Rco) and 
the total cloud radius. First, we see that the optical depths to both the H2 and CO layers 
increase with increasing radiation fields, (i.e., the transition layers are at greater column 
densities from the cloud surface). In addition, higher cloud masses result in (slightly) greater 
optical depths to the transition layers. However, Figure [7] {top) shows AAy dg to be nearly 
constant at ~ 0.6 — 0.8 over our entire parameter range. 

The cloud radii are shown in Figure [7] {bottom). The radius in CO, Rqo is fixed from 
observations by the cloud mass and column density, (eq. |2]) and is independent of the incident 



We n ote that the dens i ty in t he atomic portion of the cloud is close to the value of 230 cm ^ given by the 



theory of iKrumholz et al 



for this case (Gg = 10, Z = 1). Their result is closely tied to the present 
work, sinc e they adopt a densi ty of Srimin, where the expression for the minimum two-phase density, n„ 
is given in 



Wolfire et all (|200a) . 
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radiation field. The next largest radius, Rh2j is also found to be relatively constant with 
radiation field, reflecting the constancy of AAy oG (Fig- El top). The largest radius, -Rtot, 
which encompasses the atomic surface as well as the molecular interior, is found to increase 
slightly with increasing radiation field, a consequence of the higher column densities required 
to turn the gas molecular. We also show in Figure |8]the ratio of radii Rtot/Rco and R^^/Rco 
for cloud masses M(i?co) = 1 x 10^ M© and M{Rco) = 3 x 10^ Mq. 

In Figure M we present the calculated dark gas fraction /dg (see eq. [T]) for all of the 
cloud masses and incident radiation fields. We find the fraction is remarkably constant with 
both cloud mass and radiation field at a value of /dg ~ 0.3. As shown in equation ffTOl) . the 
constant fraction follows directly from the constant optical depth between the H2 and CO 
layers, and the constant Z' and N (or Ay) assumed in this case. 



5.2. Interpretation of Results 



The trends in Av{Rr2)j ^viRco), A^y, do, and /do can be understood from the forma- 
tion/destruction processes for H2 and CO (Appendices iB] and ICl) and the relation between 
/dg and AAv^bg (eq. [10]) • In Appendix [B| we provide an analytic treatment of Av{Rn2) 
by balancing the formation of H2 on grains with FUV photodissociation of H2. There we 
assume that the density is uniform, and since most of the mass of the cloud is in the clumps 
or, in the turbulent gas, at density (n)mcd, we use the density ric- Strictly speaking, this 
density varies somewhat from Rtot to Rn^: but since the solution heavily depends on what 
happens at Rh2, we use the value of there. The dependence on metallicity, Z' = Z/Zq, 
FUV field and density is given by equation (lB9p : 



AviRn^) - 0.142 In 



5.2 X lO^Z' 



Z'Ur 



1.75 



+ 1 



(20) 



Here, ric is in units of cm ^. We find that the column dens ity to the transition region is a func- 
tion of Gn/rir, consistent with previous investigat i ons (e.g., Sternbe rg 'l988!;_ HoUenbach fc Tielens 
I999I : IWolfire et aPboOsI : iKrumholz eraPboOsi iMcKee fc Krumholai20m 



In AppendixOwe provide an analytic treatment of y4y(i?co) by balancing the formation 
of CO by gas phase chemistry with FUV photodissociation of CO. The dependence on Gq, 
ric, and Z' is given by equation flC7|) : 



Ay(i?co) -0.1021n 



3.3 X 10' 



^0 

Z'Uc 



+ 1 



(21) 



where we have substituted the numerical values for the constants into the equation and used 
T ~ 50 K from Figure HI Taking the difference of equations fl2Tl) and fl20|) we find: 



AAv, DG = 0.53 - 0.045 In ^ - 0.097 In {Z') , 




(22) 



which leads to 



AA 



V,DG — 



0.53 -0.045 In 



) 



(23) 



for Z' = 1. Here we have assumed Gq/tIc > 0.0075Z'°-'^'^ cm'^, as is usually the case, so that 
one can ignore the factor of unity term inside the brackets of equations fl2CT]) and (^T^. As 
noted in Appendices [B] and \C\ the analytic solutions for Ay(H2) and Ay (CO) (for Z' = 1) 
agree with the numerical solutions to within 5% and 15% respectively. The analytic solution 
for AAy^DQ (for Z' = 1) agrees to within 25% of the model results (Fig. [7]). 

The principal conclusion of this analysis is that AAy ^Q is almost constant (eq. [2^ . 
consistent with our numerical results in Figure [71 Equation flTU]) shows that the dark- 
gas fraction, /do, is a function of only AAv^bg, which is nearly constant, and of Z'N or 
Ay. For conditions typical of the Galaxy — AAv,dg = 0.6 — 0.8 (Fig. [7]), Z' = 1 and 
N = 1.5 X 10^^ cm~^(or Ay — 8) — we find /do = 0.26 — 0.33, in good agreement with the 
numerical results portrayed in Figure M We now discuss the reasons behind these results. 



Why is /do so insensitive to the cloud mass? From equation (IT^ . in the limit of large 
Mach number, the median density goes as (n)med nj\4 oc na. A higher cloud mass results 
in a lower value of n oc l/M(i?co)^'^^ (eq- IHwith N22 constant), while the turbulent velocity 
scales as M{RcoY^^- Thus (n)ined oc 1/ M{RqoY/^ varies only weakly with cloud mass, with 
lower densities leading to deeper dissociation layers as the cloud mass increases. However, 
as seen in equation fl23|) . small changes in ric lead to very small changes in AAv^dg-, and thus 
to very small changes in /dg (see eq. [TOj) . 



The dark gas fraction is also quite insensitive to variation in the incident radiation field 
at constant cloud mass, M(_Rco)- For constant M{Rco), the CO radius, -Rco, and the 
distribution and value of n{r) are unchanged, although the total cloud radius increases at 
the higher radiation fields, and thus n{r) drops to lower values in the outer edges of the 



5.2.1. Insensitivity to CO Cloud Mass, M{Rco) 



5.2.2. Insensitivity to Ambient Radiation Field, Gq 



-17- 



cloud. Since AAv^bg is also relatively constant with changing Gq, the H2 radius, Rn^^ is 
unchanged. Since the radii i?H2 and -Rco are unchanged and the density distribution, n(r), 
is unchanged, we find constant masses and mass fractions. Alternatively, and perhaps more 
simply, equation f l2^ shows that AAv^dg is very weakly dependent on G'q, and therefore 
/dg is also weakly dependent as long as Ay remains fixed (eq. [10]). We note that the ratio 
of dark-gas mass to total mass (including the H I) does decrease with G'q as we add more 
shielding H I material to maintain constant M{Rqo). 

We have also investigated the effect of lowering the incident FUV and EUV/soft X-ray 
fields to G'q = CxR = 1 for the M(i?co) = 1 x 10^ Mq model. We find the dark gas fraction 
drops shghtly from /dg = 0.28 to 0.22. The H2 self-shielding becomes quite strong and 
draws the H2 transition close to the surface, resulting in Av(B.2) ~ 0.02, while ^4^(00) is 
~ 0.5. Nevertheless, the fitted functions (eqs. [20] and [5T]) behave quite well and reproduce 
the model results to within 3% and 8% respectively. 

Decreasing the incident field to G'q = 0.5, appropriate for an opaque cloud embedded in 
the ISRF, further decreases /dg to 0.18. We can check that we are getting reasonable results 
in the G'n = 0.5 c ase by comparing with observations of Nqq versus Ay- For CO dark clouds, 
Federman et al.l (jl990l ) shows CO column densities of Nqq = 4x 10^® cm~^ (twice our fiducial 



Tco = 1 surface) at total cloud columns of Ay ~ 1, and thus Ay ~ 0.5 to the CO emitting 
gas. This is reasonably consistent with our G'q = 0.5 mo del that finds Ay (R en) = 0.4. For 
H2 we note that steady state PDR models of diffuse gas (IWolfire et al.ll2008l ) find good fits 
to the observed molecular fractions deduced from UV absorption line studies. The densities 
(n ~ 30 cm~^) in diffuse gas are lower than considered here. However, applying our fitted 
formula for Ay^Rn^) to G'q = 0.5, and n = 30 cm~^, yields Ay^R^^) ~ 0.23. This is in good 
agreement with th e transition to high mo lecular fractions /(H2) > 0.1 found in the diffuse 
ISM at Ay ~ 0.26 jcillmon fc ShuliboOfih . 



5.2.3. Dependence on Cloud Column Density, N 



Equatio n f[TOl) shows t hat th e dark- gas fraction is sensitive to the mean column density 
of the cloud. iHeyer et al.l (120091 ) have argued that for Galactic mole cular clouds, the mea n 
column density within Rqq is about half the v alue we adopted from [Solomon et al.l (11987] ) , 
or N22 = 1.5/2 = 0.75. As iHeyer et al.l (j2009l ) point out, this value is quite approximate, 
since it depends on an uncertain correction for non-LTE effects in the excitation of ^^CO. 
Bolatto has found the same values for the mean column density {N22 = 0.75) and turbulent 
velocity in the SMC as Heyer adopts for the Galactic clouds. 
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To investigate the effects of variations in tlie mean column de nsity, we have also run 
N22 = 0.75 models, including the iV^/^ scaling for cr(r) suggested by lHeyer et al.l ( 120091 ) (eq. 
ITS]) . Figure ITU] (bottom) shows the low column density {N22 = 0.75) results for M(i?co) = 
1 X 10^ Mq and Z' = 0.5, 1, and 1.9. We find higher dark-gas fractions in the N22 = 0.75 
models compared to the N22 = 1.5 models; in particular, for Z' = 1, N22 = 0.75, and 
Gg = 10, the dark-gas fraction /dg ~ 0.5, or M (Ru^) / M (Rco) ~ 2. The change in /dg 
is almost entirely due to the change in mean column density; the change in a{r) accounts 
for only a few percent of the increase. There is currently no evidence for such high values 
of dark-gas mass in the local Galaxy for the hi gh mass clouds that we model here and we 
favor the higher, N22 = 1.5, value. We note that lGrenier et al.l (|2005[ ) finds higher dark mass 
fractions in local clouds with very low masses [M{Rco) < 3 X 10^ Mq]. However, these 
clouds when observed with hi gh spatial resolution CO observations revea l quite small mean 
column densities, N22 ~ 0.2 ( Mizuno et al.ll200ll : lYamamoto et al.ll2006 ) where we expect 
the dark gas fraction to be higher than /dg ~ 0.3. 



5.2.4- Dependence on Metallicity, Z' 



In Appendices |B] and O we explicitly give the dependence of the H2 and CO formation 
rates on the metallicity relative to solar, Z'. For H2, the metallicity enters due to the 
conversion from Ay to column density, N, and also in the rate of formation of H2, which 
proceeds on grain surfaces. For CO, the metallicity enters in the Ay to N conversion and 
also in the production of OH, which proceeds both on grains and in gas-phase chemical 
reactions that depend on the gas-phase abundances of oxygen and carbon. 

We ran models for our standard cloud mass, M{Rco) = 1 x 10^ Mq, incident radiation 
fields Gq = 3 ,10, and 30, and for metal abundances Z' = 0.5, 1, and 1.9. The Z' = 0.5 case 



is appropriate for molecular clouds in the Large Magellanic Clouds (LMC; lDufourlll984[ ): the 
Z' = 1.9 case is ap propriate for clouds in the Galactic Molecular Ring at a Galactocentric 
radius of 4.5 kpc (jRathborne et al" 2009 ), based on an expo nential scale length for the 
metallicity in the Galaxy of = 6.2 kpc (IWolfire et al.l 120031 ). Figure [TO] ( top) shows that 
/dg increases with lower metallicity, as predicted by equation ( ITOj) . In general, the extinction 
'^Av^BG between i^Ha and Rqo is very weakly dependent on Z' (see eq. [22]); however, the 
column density of the dark gas region goes roughly as Z'~^ so that the mass in the dark gas 
increases with lower metallicity. 

There is some indication that the dust abundance in low metallic ity systems such as the 



LMC scales as Z rather than linearly with the metallicity Z' (e.g., IWeingartner &: Draine 



2OOII ). To test the effects of this scaling we ran a model with Z' = 0.5 for gas phase metals and 
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and a dust abundance Z'^ = 0.25. We find /dg increases to /dg ~ 0.60, a value intermediate 
between a Z' = 0.25 and Z' = 0.5 scaling for botli dust and gas pliase metals (see Fig. [TOj) . 
As expected a decrease in the dust abundance at fixed N decreases Ay = 5.26Z'N22 leading 
to higher dark gas fractions; however, the increase is limited due to the metal dependence 
of Av{Rco) which leads to smaller Av{Rco) (and thus smaller AAv^dg) at Z' = 0.5 than 
at Z' = 0.25. 

We caution against extrapolating our results for varying Z' to very low values of Z'. For 
fixed N and fixed Gq, lowering Z' reduces the size of the molecular interior until first, the 
entire cloud becomes optically thin in the CO J = 1 — transition, and finally, the entire 
cloud becomes atomic H and C"*". We shall examine the dark-gas fraction for Z' <^ 0.5 clouds 
in the universe in a subsequent paper. 

5.2.5. Dependence on Visual Extinction Through the Cloud, Ay 

Equation ffTOj) shows that /dg depends only on AAv^dg/^v, where Ay = 5.26Z'iV22 is 
the mean visual extinction through the cloud. As discussed above, AAv^bg is only weakly 
dependent on Z', N221 Gq and Uc- Therefore, the main parameter that controls the dark 
gas fraction is the visual extinction through the cloud. Ay. Figure [TT] clearly shows this 
dependence; it plots /dg against Ay for two different values of Z'= 0.5 and 1.9, holding 
Gq fixed at 10 and adjusting the gas density (see eqs. HI [TH [151 aiid [18]) as N changes. 
Although, at fixed Ay, the Z' = 0.5 case does have a slightly higher value of /dg than the 
Z' = 1.9 cas^, the main parameter controlling /dg is Ay. If we fix N, lowering Z' lowers 
the Ay through the cloud, and this is the cause of the large change in /dg seen in Figure 
[TUl Essentially, lowering Ay raises the mass of the dark gas located in the surface region 
relative to the mass of the interior CO gas, which is fixed in our model. For Ay = 2, there 
is considerable shielded H2 gas out to large radius around the fixed mass CO "core". For 
Ay = 30, the CO "core" takes up much of the cloud, and the dark gas is but a thin shell on 
the surface. 

We have fixed the mass in the CO region in our models because that is what is generally 
observed. However, to better understand the dependence of /dg on Ay, it is easier to consider 
the alternate case of a cloud of fixed total mass (i.e., H I + H2 + CO mass) and variable Z'. 
As Z' is lowered, the Ay through the cloud is lowered. The radius of the CO region, -Rco, 



"'This slight change is caused by the shght dependence of AAy. dg on Z' seen in equation [22] and also 
because the different Z' cases have different N, which then means that the gas densities change in the two 
cases. 
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shrinks because of the reduced shielding in the cloud. The mass of the CO region drops, 
as Rqq since we have assumed that n oc r^^. However, the AAv^bg is relatively constant 
through the dark gas shell around the CO, so the column through the shell rises as Z'~^. 
The mass in the H2 dark gas, M(_Rh2) "~ M{Rco), changes as R^g times the column in the 
shell, where Rco < Rcs < Rn^ ■ Taking Rck = Rco to obtain the minimum mass in the dark 
gas, we see that the ratio of dark gas mass to mass in the CO region goes as Z'~^ in this 
case. Therefore, /dg increases as Ay decreases {Z' decreases) in this example. 



5.2.6. Variation with Cosmic-Rays 

The discussion of H2 and CO in the Appendices is based on the assumption that cosmic- 
ray ionization is not essential in determining the abundances of these molecules. In order to 
test the effects of higher cosmic-ray ionization rates, we start with our G'q = 10, M(i?co) = 
1 X 10^ Mq model, and enhance the cosmic-ray ionization rate by a factor of 10 from low 
columns to Ay = 2 and then drop it back to our standard value so as not to overheat the 
cloud interior. The enhanced rate tends to increase the production of He"*" from cosmic-ray 
ionization and increase the destruction of CO through reactions with He"*". The result is to 
drive the CO slightly deeper into the cloud, but to leave the H2 surface at the same depth, 
thereby slightly increasing the dark mass fraction. We find that /dg increases from ~ 0.28 
to ~ 0.39. 



5.3. Comparison With Observations 



First, we note that the dark ga s fraction fna = 0.2 8 for our standard model is in good 
agreement with the observations of (iGrenier et al.ll2005l ). who found /dg ~ 0.3. for the four 
local clouds in their sample that are more massive than M(Rco) > 3 x 10^ Mq. [f| We note 
that there is considerable scatter in their dark mass fraction for these 4 clouds, which may 
indicate scatter in the average extinctions through the clouds, but that t he average value o f 
the dark mass fraction is close to our theoretical prediction. In addition, lAbdo et al.l (I2OIOI ) 



■''Here we refer to the mass within the CO emitting regions that lGrenier et ahl (|2005l ) has designated Mhs- 
For clouds more massive than M(i?co) > 3 x 10'' Mq they found dark mass fractions, /dGj approximately 0.1 
for Cepheus-Cassiopeia-Polaris, 0.1 for Orion, 0.6 for Aquila-Ophiuchus-Libra, and 0.3 for Taurus-Perseus- 
Triangulum. For lower mass clouds they found 0.5 for Cham aeleon, 0.8 fo r Aqu ila-Sagittarius, and 0.6 for 
Pegasus. The dark gas fractions for high mass clouds found bv lAbdo et al.l (|2010l ) are 0.30 for Cepheus, and 
0.37 for Cassiopeia, and 0.36 for the low mass Polaris cloud. 
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from observations of nearby resolv e d clou ds fi nd an avera g e valu e of /dg ~ 0.34 for their 
two massive clouds. iGrenier et al.l (|2005[ ) and lAbdo et al.l (120 lOf ) find /dg higher than 0.3 
in very low mass clouds [M{Rco) < 3 X 10^ Mq] where the mean extinction is observed to 
be much less than our standard value and a higher Jdq is consistent with our prediction. 

Next, we compare our nu merical results with the observations of H I halos around 
molecular clouds in the Galaxy. IWannier et al.l (119831 ) observed H I halos around 8 molecular 
clouds and concluded they extend several parsecs beyond the CO and have a thickness of 
a few parsecs. In addition, the halos are "warm," having temperatures of at least 50 — 200 
K in order to be seen in emission over the background. We note that the use of "warm" in 
their title has been interpreted by others to mean gas at temperatures of ~ 8000 K, but in 
fact they cite temperatures of order a few 100 K. 



Andersson et al.l (Il99ll ) and lAndersson &: Wannierl (Il993l ) included additional data and 
analysis and report H I integrated intensities between 700 K km s~^ and 4300 K km s~^ 
and characteristic depths of about 4.7 pc. Based on model fits, they estimate H I densities 
nni ~ 25 — 125 cm~'^ and temperatures T ~ 50 — 200 K. We note that they suggested 
that either the formation rate of H2 was lower than the standard rate or the FUV field 
was about 10 times the interstellar value in order to simultaneously match the CO and H I 
observations. At that time they favored the lower H2 formation rate since an increased FUV 
field would produce too little OH compared to observation. However subsequent, work by 
HoUenbach et al.l (120091 ) showed that OH is formed at greater abundances than previously 
thought when grain surface chemistry is included. Thus, their suggestion that the FUV 
field is enhanced near GMCs is consistent with a similar result by Wolfire, Hollenbach, & 
McKee (2010, in preparation) in modeling the average FUV field illuminating a star-forming 
molecular cloud. 

For cloud masses in the range M{Rqo) = (1 - 10) x 10^ Mq and for Gg = 10 - 30, 
we find H I integrated intensities from 965 K km s^^ to 4100 K km s~^ and a range in H I 
halo thickness from 1 pc to 10 pc. These are in good agreement with the observations. At 
Ay[RY{^), ou r temperature range T = 70 — 80 K and density range fi = 45 — 147 cm~'^ agree 



with that of lAndersson &: Wannierl (119931 ). Therefore, we consider the model to be in very 
good agreement with these observations. 

We note that CO-to-H2 conversion factors have been inferred from ga mma-ray and far- 



infra red observations, a nd also from yirial m ass estimates. The gamma-ray (IStrong fc Mattox 



19961 ) and far- infrared (jPame et al.l I2OOII ) estimates include the "dark gas" and are X 



1.9 X 10^° cm~^ (K km s~^)~^ and X = 1.8 x 10^*^ (K km s~^)~^ from gamma-ray observa- 
tions and far-infrared observations respectively. However, these values are Galactic averages 
and observations of nearby resolved clouds show a variable X increasing towards the outer 



- 22 - 



Galaxy ( lAbdo et al.ll2010[ ). The Abdo et al. observations separate the "dark gas" and CO 
components and generally show lower X ratios [X ~ 0.87 ±0.05 x 10^*^ cm~^ (K km s~^)^^]. 
Depending on whether the dark gas is mixed with other c omponents the X val ue could be 
lower if mixed with H I or raised if mixed with CO. The ISolomon et al.l (119871 ) conversion 
factor was b ased on a virial cloud estim ate and for an idealized cloud would not include the 
"dark gas". iMcKee fc Ostriken (120071 ) re-evaluated their conversion, accounting for an 8.5 



kpc distance to the Galactic center and including the He mass, and found a revised Solomon 



et al. value oi X = 1.9 x 10^° cm ^ (K km s 



-1^-1 



In general, the X ratios derived from 



virial mass estimates exceed those bas ed on gamma-ray observations from resolved clouds 
by factors of 1.5-3.0 (lAbdo et al.ll2010l ). However, turbulent mixing of the bright CO and 
dark gas could easily cause the virial estimate to include some of the dark gas. Therefore, 
one should note that CO to H2 conversion factors are sometimes normalized to implicitly 
include the dark mass, even though it lies outside the CO gas. 



5.4. Justification of the Optically Thin Approximation 

One of the important approximations in our work is that the density fluctuations in 
the molecular cloud — the clumps — are optically thin, or Ay < 1 through each clump. This 
approximation can now be justified in light of our results: the fact that the dark-gas fraction 
depends primarily on Ay, the mean visual extinction through the entire cloud, but only 
weakly on the density and radiation field, implies that it is insensitive to the distribution of 
matter within the cloud. Consider an extreme example of a cloud with very opaque clumps, 
namely a cloud with clumps so opaque that their mean extinction is equal to that of the 
entire cloud. Ay ~ 8. We assume that the cloud with opaque clumps has the same mean 
Ay, and clump density, ric, as the cloud with transparent clumps. For simplicity we assume 
that there is no radial variation in the properties of the clumps. Let a be the clump radius. 
Then the number of clumps. A/", is determined by the condition that the clumps occupy a 
fraction fy of the volume: 

Afa' = fyRl,. (24) 
The condition that the clumps have the same Ay as the cloud as a whole gives 

rica = nRtot = /v^c-Rtot, (25) 

so that a = /y-Rtot- When viewed from the outside, the fraction of the cloud that is covered 
by clumps (the projected covering factor) is then 



(26) 
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as expected. From the perspective of a clump inside the cloud, this means that about half 
the sky is covered by clumps, since, on average, the path length from a clump to the surface 
of the cloud is half the path length through the entire cloud. The average radiation field 
incident on a clump is therefore reduced by about a factor 2 compared to the radiation field 
incident on the cloud whole. 

We can now estimate the dark-gas fraction of a cloud of these opaque clumps. Each 
clump will have a dark-gas fraction that is almost the same as that of the entire cloud, 
modified only by the reduction in Gq by about a factor of 2. Since the extinction of the dark 
gas layer depends only weakly on the intensity of the radiation field, AAv^bg oc 0.045 In Gq, 
this modification is very slight for typical values of Ay ~ 8 (see eg. [TOj) . This argument does 
not prove that clumps with Ay intermediate between the large value ~ 8 assumed here and 
the small values Ay < 1 assumed in the rest of this work will also have the same dark-gas 
fraction. However, because this range 1 < Ay < 8 is small it strongly suggests that this is 
the case. We conclude that inclusion of finite optical depth of clumps in a cloud is unlikely 
to significantly alter our conclusions. 



5.5. Justification of Assumption of Steady State Chemistry 

Our models solve for the steady state abundances of the atomic and molecular species. 
In this subsection we examine the assumption of steady state chemistry in the turbulent 
molecular cloud surfaces by comparing the chemical and dynamical times at the cloud depth 
A{Rii^). The chemical timescale, tchem, is the time for atomic gas to become completely 
molecular {rin^ = 0.5(n)med), thus tchem = ^-5/ (TZ{n)med) where 71 = 3 x lO'^'^Z' cm^ 
is the rate coefficient for H2 formation on dust grains (see Appendix [B]) . To compute the 
dynamical timescale, tdyn, we require the characteristic distance for the turbulence to bring 
molecular gas from the interior to the surface, and to bring atomic gas from the surface to 
the interior, where the interior is Ay greater than AylRn^)- This characteristic distance 
is ddyn = 1-9 X 10"^^ Ay {Rii2)/n. The dynamical timescale is this distance divided by the 
turbulent velocity for lengthscale (idyn, ^dyn = '^dyn/o"(c?dyn) where a is the ID turbulent 
velocity dispersion (eq. [T8|) . 

For a cloud mass M{Rco) = 10^ Mq, and Z' = 1, we find n ^ 140 cm^^, and (n) ^ 
600 cm"'^ at Ay{Rii^) = 0.22 and thus tchcm ~ 2.8 x 10^^ s, c^dyn ~ 1.0 pc, and tdyn ~ 
4.3 X 10^'^ s. Comparing chemical and dynamical timescales we find tchem/^dyn = 0.65. 
Similarly, for M(i?co) = 10^ Mq we find tchem/^dyn = 0.48. For tchem/^dyn ~ 1, a steady state 
approximation is valid since gas has time to reach steady state before significant turbulent 
mixing. We find ratios less than one, but only marginally so, and thus we expect some affects 
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from turbulence. We note however that turbulence brings atomic gas inward past Av{Rr2) 
into the molecular gas, and at the same time it brings molecular gas to the surface. Thus, 
the first order effect is to spread out the rise of xhj, but not to move Av{Rn2)- addition, 
we have demonstrated in subsection (15. 3p that the steady state model is in good agreement 
with observations, and thus the effects of turbulence are modest. 



6. Discussion and Summary 
6.1. Model Assumptions 



We have constructed models of molecular clouds to investigate the fraction of gas that is 
mainly H2 and contains little CO. These conditions exists where the CO is photodissociated 
into C and C"^ but the gas is molecular H2 due to either H2 self-shielding or dust shielding. 
Such conditions can exist either on the surfaces of molecular clouds or the surfaces of clumps 
contained within such clouds. Observations indicate that the mass in this "dark gas" can 
be as high as 30% of the total molecular mass in the local Galaxy. Previous theoretical 
plane-parallel models of ind ividual PDRs have indicated that such a layer should exist (e.g., 
van Dishoeck fc Blacklll988l ). but here we construct models directed towards molecular clouds 
as a whole while including observational constraints on cloud mass, radius, average density, 
and line width and theoretical considerations of the likely strengths of the UV fields impinging 
on GMCs. 

We assume that the surface of each cloud is isotropically illuminated over 27r steradians 
by a soft X-ray/EUV field and an FUV radiation field. We use the standard cosmic-ray 
ionization rate of 1.8 x 10~^^ s^^ per hydrogen nucleus everywhere in the cloud for all cases 
but one test case. There is evidence from observations of Ht that the cosmic-ray ionization 



rate is a factor of 10 higher in some portions of the diffuse ISM (llndriolo et al.ll2007h. bu t 
there is no indication that such rates apply in molecular cloud interiors (IMcCall et al.lll999l ). 
In preliminary work of M. Wolfire et al. (2010, in preparation) we find that the average FUV 
field on clouds is ~ 20 times the local Galactic interstellar field. This elevated field arises 
from the distribution of OB associations around the cloud. 

The distribution of temperature, densit y, and abundances wit hin the HI, H2, and CO 
layers are calculated using the FDR code of iKaufman et al.l (120061). In constructing model 
clouds from the FDR output, we impose the constraint from ISolomon et al.l (119871 ) that 
the typical column density through the cloud is N22 = 1.5, independent of cloud mass and 
independent of radius inside a given cloud (since n oc r~^). The locally averaged density n 
and radius of the CO cloud, -Rco, as functions of the CO-cloud mass, M(/2co) follow from 
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equations ([2]) and (jl]). In our notation, M(i?co) is the molecular mass contained within the 
tqo = 1 surface (of the CO J = 1 — transition) including the mass of H2 and He, and Rqo 
is the radius of this CO photosphere. We impose a volume-averaged density distribution 
n{r) that behaves as n oc 1/r throughout all layers of the cloud. Note that in regions where 
the thermal pressure Pth lies between Pmin and Pmax both warm (T > 7000 K) and cold 
(T ^ 500 K) gas solutions are possible. In this regime we use the cold solution from our 
model results with ric the density of the cold clumps and n the average over cold and warm 
gas components. The warm gas fills the volume, but contains little of the mass. 

We test two extreme limits for clumps within the cloud, one in which all clumps are 
optically thin to the incident FUV field, and one in which all clumps are optically thick to 
the incident FUV field. In the optically thin approximation, the PDR model output as a 
function of Ay gives directly the distribution in Ay throughout the cloud. In the optically 
thin limit we test two different models for the cloud density distribution. First we use 
constant thermal pressure models and second we include two sources of thermal pressure, 
radiative heating (which by itself would lead to a two-phase equilibrium) ar id supersonic 



turbulence. The distribution of two-phase thermal pressure is calculated as in lWolfire et al. 



(120031 ) and stored in a look-up table as a function of total column density and molecular 
fraction (see Fig. [2]). This pressure drops as one moves into the cloud due to the absorption 
of the radiation responsible for heating the gas. Supersonic turbulence is char acterized by a 



mass -weighted median density {n)^ed = nexp{fi), with fi = 0.5 ln(l-|-0.25A^^) (jPadoan et al. 
19971 ). When the two-phase pressure drops below P^^'^^ = Xt{n)^Q^kT , where Xt is the sum 



over the abundances of all species relative to hydrogen nuclei, we assume that the turbulence 
maintains the gas at a thermal pressure P^^^^- The sound speed that enters in the Mach 
number is calculated from the PDR model output while the turbulent velocity is given by 
the linewidth-size relation (eq. [H 



In the limit of very optically thick clumps (with optical depths at least as large as the 
average optical depth through the cloud), we find that the average radiation field on a clump 
ranges between about G'q/2 to G'q. In this limit, the dark mass fraction of a clump is the 
same as the dark mass fraction of the entire cloud. Since the dark gas mass fraction of a 
spherical clump or a spherical cloud is insensitive to the incident FUV field (see eqs. [TD] and 
1221) . the dark mass fraction of the cloud does not change significantly compared with the 
case of a cloud made up of optically thin clumps. 

We have tested the steady state assumption for the chemistry by comparing the time 
to form molecular hydrogen, tchem, with the dynamical time, tdyn, for turbulence to bring 
molecular gas from the interior to the surface and to bring atomic gas to the interior. We 
find tchem/^dyn ~ 0.7 for a cloud mass of M(i?co) = 1 x 10^ M© and tchemAdyn ~ 0.5 for a 
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cloud mass of M(i?co) = 1 x 10^ Mq. Thus we expect modest affects due to turbulence; 
mainly to spread out the transition from atomic to molecular hydrogen but not to move 
Av{R}i2)- We also note that steady state models agree well with observations. 



6.2. Dark- Gas Fraction 



For our standard cloud mass, M{Rco) = 1 x 10^ Mq, and incident radiation field, 
Gq = CxR = 10) W6 find dark-gas mass fractions of /dg = 0.28 and /dg = 0.31 for constant 
thermal pressures of Pth = 10^ K cm~^ and Pth = 10^ K cm~^, respectively. These correspond 
to a total molecular mass M(i?H2) ~ l-4M(i?co)- For models that include both thermal 
and turbulent pressures, we find essentially the same results as for the cases with constant 
thermal pressure. The variation in /dg ranges from 0.25 to 0.33 over a range in Gq from 
3 to 30 and a range in cloud mass from 10^ Mq to 3 x 10^ Mq (Fig. [9]). Figure |3] shows 
the distribution in thermal pressures and densities and Figure S] shows the distribution in 
temperatures and chemical abundances for the standard model. 

The constant value of /dg for fixed Ay can be understood from the analytic solutions for 
Av{Rk2) S'lid Av{Rco) in Appendices [B] and [Cl and the expression for /dg in equation ffTOj) . 
In the limit of G'^/n > O.OOTSZ'"-^^ cm^, both AyiRn^) and Av{Rco) increase as ln(G"o/ra). 
Thus, the optical depth through the H2 layer AA^^dg is a weak function of G^/n and Z' 
(eq. [22]) and is nearly constant over our parameter space. Furthermore, we find that /dg 
is a function of only AAv^bg and the mean extinction through the cloud Ay = 5.26Z'N22 
(Eq. [TOj) : thus, for a given Ay, the dark-gas fraction is constant. However, /dg increases 
significantly if Ay decreases. 

Our nume rical results compar e well with observations of the local Galactic dark-gas 
fraction, which iGrenier et al.l (120051 ) find to be /dg ~ 0.3, when averaged over their four most 
massive clouds with masses between 3 x 10^ - 3 x 10^ Mq. Lower mass clouds observed by 
Grenier et al.l (|2005[ ) are observed to have low Ay and thus high dark-gas fractions consistent 
with our prediction. Our H I integrated intensities of 965 K km s~^ to 4100 s~^ , H I cloud- 
halo thickness of 1 pc to 10 pc, average cloud densities of n ~ 45 — 150 cm~^, and average 
H I temperatures of ~ 70 — 80 K are in good agreement with t he ob servations of H I cloud 



halos observed by IWannier et al.l (1l9831 ) , lAndersson et al.l (Il99ll ) , and lAndersson fc Wannier 
f|l993h . 



We have carried out several additional tests to assess the dependence of our results on 
the cosmic-ray ionization rate, the clump optical depth, the mean cloud column density, the 
metallicity, and the mean visual extinction through the cloud. We ran our standard model 
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with the cosmic-ray ionization rate a factor of 10 higher at Ay < 2, as suggested by Hjj' 
observations in some regions of the diffuse ISM, and found that this enhanced cosmic-ray 
rate only shghtly increases the dark-gas fraction. In the limit of very optically thick clumps 
we find no significant change in /dg- 



Heyer et al.l (120091 ) has suggested that the mean column density (and the refore the mean 



visual extinction Ay) in Galactic CO clouds is about half the value found by ISolomon et al. 



(jl987l ). i.e., N22 — 0.75 or Ay — 4. The results of changing the mean column density, the 
metallicity, and the mean visual extinction through the cloud are illustrated in Figures [TD] 
and [m We find that the dark-gas fraction increases with lower extinction since the dark 
gas occupies a larger fraction of the cloud (see eq. [10]). There is also a weak dependence 
of AAy^DQ on the mean column density that tends to slightly mitigate the dominant effect. 
Lower mean columns lead to lower mean densities (eq. H]) and lower ric (eqs. and [T^ . and 
thus lower AAy^^c (eq- [22]). We examine metalhcities appropriate for the LMC {Z' = 0.5), 
for the local Galaxy {Z' = 1), and for the molecular ring at i? = 4.5 kpc {Z' = 1.9). In 
general, /dg increases as the metallicity drops for fixed columns because the mean extinction 
through the cloud decreases, which raises the ratio of the surface dark gas to the interior CO 
gas (again, see eq. [TU]). There is also a weak dependence of AAy ^Q on Z' that also slightly 
increases AAy^^c (or /dg) with decreasing Z' (eq. [22]) even if the column is changed so that 
Ay remains fixed. We note that in the case of varying N and Z', but at constant NZ' or 
Ay (see Figure [11]), the change in /dg is entirely due to the weak dependencies of AAy dg 
on N and Z' as noted above and shown in equation f[221) . We also examine the case for dust 
scahng as Z''^ (= 0.25) while gas phase metals scale as Z' (= 0.5). We find /dg is larger 
than when both metals and gas scale together. 

In Appendices [B] and [C] we derive analytic solutions for Ay{Rii^) and Ay{Rco) as 
functions of density n, FUV field Gq, and metallicity Z'. For the case of H2 we find an 
expression for the abundance of H2 by balancing formation and destruction processes and 

0.25n. Similar s tudies have been carried out by. 
We find our fits are good to 



Sternberg 


( 


1988 


); 



±5% for all models presented in this paper. The fits to CO are found by integrating the 
expression for the abundance of CO to a column density of A^(CO) = 2 x 10^^ cm~^, where 
Tco = 1 (the optical depth of the CO J= 1-0 transition). The fits are generally good to 
within ±15% except for the lowest Z' and Gq model, where the fit is good to ±25%. We 
have also derived analytic expressions for AAy^ dg and /dg in the main text (eqs. [22] and 
[TUl respectively). 

The overall result of this paper is a theoretical derivation that /dg — 0.3 ± 0.08 for 
GMCs with Ay ~ 8 in our Galaxy. Therefore, a significant fraction of the molecular gas 
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in our Galaxy lies outside the CO gas. As discussed in §2.4.2, some calibrations of the 
CO line intensities to H2 mass take into account this "dark" H2 gas. However, it is im- 
portant to be aware of this component since it contributes significantly to the gamma ray, 
infrared/submillimeter continuum, and [CII] 158 /xm emission from clouds in galaxies. Its 
contribution to the star formation in a galaxy is as yet undetermined. The importance of 
this component increases as the metallicity decreases, such as in the outer regions of galaxies 
or in low metallicity galaxies. 

Although the gas in the C+/H2 layer is termed "dark gas," it emits [C II] 158 /im line 
emission and the dust in the dark layer emits infrared continuum. In a subsequent paper we 
will estimate the emission from the "dark gas". In addition, a future paper will model in 
detail the lower metallicity clouds found, for example, in the SMC and early universe. 
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A. Dark-Gas Mass Fraction 



Here we determine the dark-gas mass fraction, 

M{Rn,) - M{Rco) 



/] 



DG 



M(i?Hj 

for clouds with a power-law gradient for the mean density. 



(Al) 



(A2) 



for kp < 3 (the reason for no rmalizing the density at R^^ will become apparent below). 



GMCs typically have fcp ~ 1 ( lLarsonlll98ll ). which is the case considered in the text. For 
kp < 3, the mass inside a radius r is 



M(r) 



3 kf, 



and the mean column density inside r is 

M(i?Hj 4:RnSiRn,) 



N{Ru, 



Ha 



3 kn 
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where /xh is the mean mass per H nucleus. For kp 1, the column density of the dark-gas 
layer is 
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(A5) 
(A6) 



where we used equation ( lA4p in the second step. The mass fraction of the dark gas (eq. lAll) 
becomes 

'^tco 
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with the aid of equation (1A3[) . Equation (1A6P then implies 
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(A9) 



3-kJ Z'N22{Rn,] 

where we used equation ([3]) in the second step. Taking the limit as fcp — )■ 1 gives equation 
(IH]) in the text, since AA^dg = NqAAv^dg and N is constant, so the argument i^Ha is not 
needed. 

In order to see how the dark-gas fraction depends on the density gradient, we rewrite 
equation (1A9P as 
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If the dark gas corresponds to a relatively small fraction of the total extinction [x <^ 1), 
then 



X 



/dg(M - fl)G{kp = 1) - ^• 



(A13) 



For < A;p < |, where \a\ > 3, the magnitude of the difference between the actual dark- 
gas fraction and that found in the text, where kp is taken to be unity, is < This is 
typically small: For AAv^dg < 1 and Z'N22 — 1-5, as in the text, we have x < 0.5 and 
|/dg(^p) - /dg(^p = 1)1 < 0.04. For the case N{Rh2) < 4AA^dg {x > 1) we see from 
equation flAlOP that /dg ~ 1- 
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B. Analytic Treatment of Av{B.2) 

The analytic treatment of the H/H2 transition assumes that H2 is formed on interstellar 
dust grains with an effective formation rate coefficient of 7?. = 3 x 10~^^ cm^ s~^, and is 
destroyed by FUV photodissociation with an unshielded rate of Jq = 6 x 1 0~^^ s~^ in the 



Habind (119681 ) field (1.6 x lO^^ erg cm'^ s'^) or = 1.02 x 10-^° s'^ in the iDraind f|l978l ) 
field, which is 1.7 times stronger. In order to simulate an isotropic field incident from 27r 
steradians on an optically thick slab, we assume a ID FUV flux incident at an angle of 60 
degrees with the normal to the slab. Consequently, to penetrate to a (normal) distance x 
into the slab, a photon needs to travel 2x. In a steady state in which the photodissociation 
of molecular hydrogen is balanced by the formation on dust grains, we have 

G'o/^/,nH,e-2^«2^- = nZ'nnni, (Bl) 



where G'q is the energy density of the dissociating radiation fleld in units of the lDraind (119781 ) 
fleld, fs is the H2 self-shielding factor (see below), n-n^ is the number density of molecular 
hydrogen, b^^ is the dust attenuation factor that accounts for differences in attenuation be- 
tween the visible and the FUV and also approximately includes the effects of FUV scattering. 
Ay is the visual extinction into the slab normal to the surface (see eq. E]), is the hydrogen 
nucleus number density, and tt-hi is the atomic hydrogen number density. In this Appendix 
and the following one, we are treating the gas in the clumps, where most of the mass re- 
sides. Therefore, there is no need to distinguish the density n from the local clump value Uc- 
The factor of 2 in the exponential follows from the factor of 2 longer pathlength due to the 
oblique incidence of the FUV flux (see above). We assume that the surface area of grains 
increases proportionally to the increase in metallicity Z', so that both the H2 formation rate 
coefficient and the FUV extinction increase linearly with Z'. 

In order to obtain a relatively simple analytic approximation to the variation of riHa with 
depth into the slab, we require that the self-shielding factor be approximated by a power 
law with the form 



fs 



(B2) 



_2iVH,. 

where Ni is a constant that is determined by comparison with detailed numerical calculations 
(see below), A'h2 is the H2 column density measured normal to the cloud surface, and the 
factor 2 is due to the longer pathlength associated with the 60 degree incidence angle. 

We wish to flnd an analytic expression for Ay{Rii2), where Ay^Rn^) is deflned as the 
characteristic value of the (normal) visual extinction to the point in the slab at which the 
hydrogen makes the transition from atomic to molecular form. We deflne this to be at the 
point where riu^ = 0.25n and riui = 0.5n. In other words, Ayi^Rn^) is where the hydrogen is 
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half atomic and half molecular]^ We denote the corresponding column densities by N{Rii^); 
for example, Nn^^Rn^) is the column of H2 from the surface of the cloud at Rtot to Rh2, 
the column of H2 to a depth given by Av{Rh2)- 

Substituting riHa = 0.25n and nni = O.bn into equation (IBip we obtain for the H2 
column A^H2(-Rh2) corresponding to Av{Rr2) 

We can also multiply equation (1B1|) by dx and integrate into the slab to obtain an 
expression for A'^h2 ^ls a function of or Ay into the slab. We obtain the general solution up 
to the point A^H2(-Rh2) or AyiR^^) where the gas is half atomic and half molecular. In order 
to simplify this integration, we assume that whi = c\n and is constant for < A^H2(-Rh2)- 
Note that C\ varies from 1 at small column to 0.5 as we approach Ay[RYi-^^ and it can be 
adjusted to best match the numerical code results. Since the solution heavily depends on 
what happens near AyiR^^, we expect c\ ~ 0.5; in addition, we expect the density to be 
n ~ n(i?H2)- For d <\ (which we show later is always the case) we obtain: 



A^H^ = 0.5 



d 



Substitution of Ay = Ay{Ru2) into this equation yields another expression for iVH2(-RH2) as 
a function of Ay{R}i2)- 

Equating the two expressions for A^H2(-Rh2) (eqs. lB3l and IB4p . we obtain a trancendental 
equation whose solution gives Ay{Rii^): 



- d)cinZ'nNo 



■2bYi2Av{RYi2)/d 



(B5) 

We can obtain analytic solutions in the two limits of small Ay(i?H2) and large Ay{R}:{^). In 
the hmit Ay^R^i.^) <^ [Zh-^^Y^ ~ 0.25 (see below), which corresponds to no appreciable dust 
extinction of FUV and entirely self-shielding by H2, we obtain 



®Note that the simple power law expression given in equation ()B2p need only apply for shielding columns 
Nyi2 ^ 10^^^^° cm~^ that correspond to typical H2 columns at positions in the slab with Ay somewhat less 
than or equal to Av{R]i2)- 
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In the limit Av{Rh2) ^ [2&H2] ^ ~ 0.25, which corresponds to significant dust extinction as 
well as self-shielding, we obtain 



Ay(i?H2) 
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1 - d)ciNQ yinz'n 
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(B7) 



These limiting solutions correspond to Gg/n -C 0.019 cm^ and (?o/n ^ 0.019 cm^, respec- 
tively. The two limiting solutions can be combined in a simple expression that preserves 
the limiting solutions, and smoothly transitions from one to the other when Av{R]i2) ~ 
[26h2]"^ ~ 0.25: 
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Our numerical code uses the Meudon code to solve for the dust extinction and H2 self- 
shielding that leads to the H/H2 transition in a slab. From the code results for Z' = 1, we find 
that for our range of G'^/n (which extends from 4.2x 10~^ cm^ to 0.3 cm^), Ay{R}i,^) ~ 0.05—1 
and iVH2(-RH2) ~ 10^^ — 10^° cm"^. We fit the self shielding for the critical range 10^^ 
cm-2 < < 10^° cm-2 with A^i = 3.6 x 10^^ cm-^ and d = 0.57. We fit the dust shielding 
of the FUV with 6h2 — 2. In the hmited allowed range 0.5 < Ci < 1, we find that Ci = 0.5 
best fits the numerical results, as expected from the above discussion. Therefore, we obtain 
the analytic expression 



AviRn^) ^ 0.142 In 



1 + 5.2 X lO^Z' 



Z'n 



1.75 



(B9) 



where n is in cm~^ and references the Drainc field. Over the relevant range 4.2 x 10~^ < 
Gq/ti < 0.3 (which corresponds to 0.05 < Ay{RYi,^) < 1.1), and for Z' = 0.5, 1, and 1.9, we 
find that this analytic solution matches the numerical code result to better than 5%. 



C. Analytic Treatment of yly (CO) 

In our standard runs with relatively low cosmic-ray ionization rates (primary rates 
~ 2 X 10~^^ s~^), the chemistry of CO from the surface of the cloud to ^4^(00) is dominated 
by the formation of OH on grain surfaces followed by a chemical chain that leads to CO. 
The OH abundance is determined by equating the formation of OH on grain surfaces (O -|- 
gr ^ OH) to the FUV photodissociation of OH (OH + hiy ^ O + H): 

Tixon^ = lUG',e-''°^^^xoKn, (CI) 
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where n is the hydrogen nucleus density, 71 ~ 5 x 10~^^Z' cm^ s~^ (IHoUenbach et al.ll2009[ ) 
is the effective rate coefficient for O atoms coUiding and sticking to dust grains (the dust 
cross section per H nucleus is contained in 7i), xnn = n(0) is the den sity of O atoms in 



the gas, = 3-5 x 10-^° s'^ flRoberge et al.lll99ll : IWoodall et al.ll2007l ) is the (unshielded) 



photodissociation rate of OH in the local (Draine) interstellar field, Gq is the ratio of the 
incident FUV field to the Draine field, 6oh = 1-7 is the factor that, multiplied by Ay, gives the 
effective grain optical depth in the FUV to photons that dissociate OH, xqh^ = n(OH) is the 
OH number density in the gas, and the factor of 2 in the exponential reflects the assumption 
that a diffuse field incident on a cloud can be approximated by a ID flux incident at an angle 
of 60 degrees to the normal. Note that we assume that every atomic O that strikes a grain 
sticks to the grain, reacts with an H atom, and comes off the grain as OH. Equivalently, it 
might react with another H atom, form H2O on the grain surface and photodesorb as either 
OH or as H2O, which immediately photodissociates to OH. We derive from eq. fICip : 



XOK 



ro^G'oe-^''onAv 



(C2) 



Once OH is formed the chemistry proceeds as follows: 
OH + C+ ^ C0+ + H. 
C0+ + H2 ^ HCO+ + H. 
HCO+ + e ^ CO + H. 



CO + hu ^ C 



O. 



The first reaction dominates the formation of CO^; the second reaction dominates the de- 
struction of CO"*" and the formation of HCO'*'; the third reaction dominates the destruction 
of HCO"*" and the formation of CO; the last reaction dominates the destruction of CO. As a 
result, every C0+ that is formed by the first reaction results in a formation of a CO by the 
third reaction. We can then equate the formation rate of CO (the first reaction) with the 
destruction rate of CO (the last reaction): 



72a:oHa:c+^ 



/coG"o/coe~''^°^-xcori, 



(C3) 



where 72 = 2.9 x 10^9(7/300 K)"^-^^ cm^ (Dubernet et al. 1992; E. Herbst 2006 private 
communicatioi]0) is the rate coefficient for OH reacting with C"*" (first rea ction), xc+n is the 



C"*" gas phase number density, I'^q = 2.6 x 10 s ^ ( IVisser et al.l 120091 ) is the unshielded 



photodissociation rate of CO in the Draine field, fco is the shielding factor caused by CO 



^see also http://www.physics.ohio-state.edu/~eric/research.html 
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(self-shielding) and H2, the exponential factor is shielding of CO by dust, 6co = 3.2, and 
xcon is the number density of CO. 



Near Av{Rco), the CO photosphere, we find by fitting to IVisser et al.l (120091 ) that 

-0.60 



fco — Ceo 



2N, 



CO 



(C4) 



where cqo = 4.4 x 10~^ and where the factor of 2 applies because A'^co is the column of 
CO along a line of sight perpendicular to the surface, whereas our isotropic field is roughly 
equivalent to a ID flux incident at an angle of 60 degrees to the surface. The flux then 
traverses a column 2Nco to the point in question (see the discussion of Av{Ru2) Appendix 



We substitute equation (IC4p and equation (lC2p into equation (1C3P and integrate both 
sides of the equation over dz, the perpendicular depth into the cloud, until we reach the 



CO photosphere, N( 



CO = 2 X 10^6 ^^-2 AviRco) 
dNco: 



Here we use ndz = N^dAv with 



A^o = 1-9 X 10^7^' and xcondz 



-'oh'-'o 



■2x1016 cm-2 



Tl pi 

-'co'-^o 



2A^< 



-\ -0.60 



CO 



1016 



cm 



dN< 



CO- 



(C5) 
(C6) 



Assuming that the density is constant from the surface to y4y(i?co) and that the carbon is 
entirely C"*" (xc+ = xq), we solve for Ay (CO): 



AviRco) 



1 



2(&oH + &co' 



In 



12 



4.35 X 10i6cco(6oH + &co)/oH-^c oG^o 
A'o7i72a;oa;c«^ 



+ 1 



(C7) 



where xq ~ 3.2 x 10~'^Z' is the gas phase abundance of atomic O, xq — (1.6 x 10~^)Z' is the 
gas phase abundance of atomic carbon, and Z' is the metallicity relative to solar. We find 
this analytic solution is good to within 15% except for the lowest Z' and G'q model {Z' = 0.5, 
Gq = 3), where the fit is good to within 25%. 
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Figures 




Fig. 1. — Illustration of a model cloud showing the radius it!co of the CO core, the radius i^Ha 
where 2nn2 — '^hi (equal mass density in H atoms and H2 molecules) , and Rtot the total cloud 
radius. Within R < Rco, gas is mainly CO and H2. Within the range _Rco < R < Rn2, 
gas is mainly H2 whereas the gas phase carbon is mainly C and C^. Within the range 
i?H2 < R < Rtot gas is mainly H I whereas the gas phase carbon is mainly C^. Ay{R}i,^) is 
the optical depth measured from the outer radius to Rn^ and Ay (Rco) is the optical depth 
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measured from the outer radius to i?co- 
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Fig. 2. — Average two-phase thermal pressure P^^ = {PmmPmax)^^^ as a function of total 
column density from the cloud surface N and molecular fraction /(H2) = 2Ny{^/N . The FUV 
radiation field and soft X-ray/EUV ionization rates are a factor of 10 times higher than local 
interstellar medium values. Curves are shown for /(H2) = {solid curve), /(H2) = 0.1 
{dotted curve), and /(H2) = 0.5 {dashed curve) 
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Fig. 3. — Top panel: Thermal pressure as a function of optical depth into the cloud for 

M{Rco) = 1 X 10^ Mq, N22 = 1.5, Z' = 1, and incident radiation field G'q = Cxr = 10- 
Curves are shown for two-phase thermal pressure P^^/k = Xtn2pT/k {dotted curve) and for 
turbulent thermal pressure P^^^y^ = Xt{n)„^cdT/k {solid curve). Optical depths at Ay(i?H2) 
[nHa/''^ = 0.25] and Av{Rco) (tco = 1) are indicated by vertical dashed lines. Thermal 
pressure in the turbulent medium dominates at Ay > 0.001. Bottom panel: Density as 
a function of optical depth into the cloud for M{Rco) = 1 x 10^ Mq, N22 = 1.5, Z' = 
1, and incident radiation field Gq = Cxr — 10- Curves are shown for two-phase density 
n2p, mass- weighted median density (n)nied in a turbulent density distribution, and volume- 
averaged density n. The two-phase density distribution is shown only where two-phase 
pressure dominates {Ay < 0.001). The local (model) density ric = n2p where two-phase 
pressure dominates and = {n)^cd where turbulent pressure dominates. Optical depths at 
Av{Rh2) [^n^/^ = 0.25] and Av{Rco) {^co = 1) are indicated by vertical dashed lines. 
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Fig. 4. — Top panel: Temperature as a function of optical depth into the cloud for M{Rco) — 
1 X 10^ Mq, N22 = 1-5, Z' = 1, and incident radiation field Gq = Cxr = 10- Bottom panel: 
Abundances of C"*" {long-dash curve), C° {solid curve), OH {dotted curve), and CO {short- 
dash curve) as functions of optical depth into the cloud for M{Rqo) = 1 x 10® Mq, N22 = 1-5, 
Z' = 1', and incident radiation field Gq = Cxr = 10- Note that we do not include freeze out 
of H2O on grain surfaces, which would affect OH abundances at Ay > 3. 
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Fig. 5. — Top panel: Thermal pressure as a function of optical depth into the cloud for 

M{Rco) = 1 X 10^ Mq, N22 = 1.5, Z' = 1, and incident radiation field G'q = Cxr = 10- 
Curves are shown for two-phase thermal pressure P^^/k = Xtn2pT/k {dotted curve) and for 
turbulent thermal pressure P^^^y^ = Xt{n)„^cdT/k {solid curve). Optical depths at Av{Rh2) 
[n-a^/n = 0.25] and Av{Rco) (tco = 1) are indicated by vertical dashed lines. Thermal 
pressure in the turbulent medium dominates at all Ay- Bottom panel: Density as a function 
of optical depth into the cloud for M{Rco) = 1 x 10^ Mq, N22 = 1.5, Z' = 1, and incident 
radiation field Gq — (^r — 10- Curves are shown for mass-weighted median density (n)i„ed 
in a turbulent density distribution, and volume- averaged density n. Turbulent pressure 
dominates at all Ay and thus no two-phase density is shown. The local (model) density 
nc = {n)mcd- Optical depths at Av{Ru2) [^^2/1^ — 0.25] and Av{Rco) (tco = 1) are 
indicated by vertical dashed lines. 



-47- 



1000 



100 



E- 



10 



1 1 1 1 1 M 1 1 1 1 1 1 1 1 M 1 1 1 1 1 1 1 M 1 1 


1 1 1 1 1 M 1 1 1 1 1 1 1 1 L 




: U{RJ = IxlO^ M,, 








1 1 1 1 


1 



1 

0.0001 0.001 



0.01 0.1 



10 



0.001 
0.0001 




10-11 



A,(R„ ) A,(R,,) 



_i I I I I I I I I I I I I I 



0.0001 0.001 



0.01 0.1 

A. 



10 



-48 - 



Fig. 6. — Top panel: Temperature as a function of optical depth into the cloud for M{Rco) — 
1 X 10^ Mq, N22 = 1-5, Z' = 1, and incident radiation field Gq = Cxr = 10- Bottom panel: 
Abundances of C"*" {long-dash curve), C° {solid curve), OH {dotted curve), and CO {short- 
dash curve) as functions of optical depth into the cloud for M{Rqo) = 1 x 10^ Mq, N22 = 1-5, 
and incident radiation field G'q = Cxr = 10- Note that we do not include freeze out of H2O 
on grain surfaces, which would affect OH abundances at Ay > 3. 



-49 - 



1 1 


1 1 1 1 1 1 

M(Reo) = 


1 1 1 

3x106 Mq 


1 1 1 1 1 

(top) J 


; A,(R„^) 


M(Rco) = 


1x105 Mq 


(bottom) J 


: ^^v.DG 








1 1 — 


; rrT\ 1 


1 1 1 






1 10 100 



1 1 1 


1 1 1 1 1 1 1 

M(Rco) 


1 1 1 

= 3x106 


1 III- 

Mo 










! ^lot 










M(Rco) 


= 1x105 


Mo - 










1 1 1 


1 


1 1 1 


1 1 1 1 



1 10 100 

Go' 



-50- 



Fig. 7. — Top panel: Optical depth as a function of cloud mass, M(i?co), and incident 

radiation field, Gq = Gq/1.1. Curves are shown for N22 = 1-5, Z' = 1, and two cloud masses, 
M{Rco) = 3 X 10^ Mq and M{Rco) = 1 x 10^ Mq. The optical depth from the cloud surface 
Rtot to Rco is ^y(-Rco) {dotted curve) and the optical depth from the cloud surface i?tot 
to i?H2 is Av{Ru2) {solid curve). Also shown is AAv,dg = Av{Rco) — Av{Rr2) {dash-dot 
curve) Bottom panel: Cloud radii i?tot {dashed curve), Ru^ {solid curve), and Rco {dotted 
curve) are shown as functions of cloud mass M{Rqo) and incident radiation field normalized 
to the local interstellar field, Gq = Gq/1.1. Curves are shown for N22 — 1-5, Z' — 1, and two 
cloud masses M{Rco) = 3 x 10^ Mq and M{Rco) = 1 x 10^ Mq. The radius Rn^ is where 
half the nuclei are in H2 {n-n^/nc — 0.25), and the radius Rco is where rco = 1- 
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Fig. 8. — Ratio of cloud radii Rxot/Rco {dashed curve), and R^i^/Rco {solid curve) are 
shown as functions of cloud mass M{Rqo) and incident radiation field normalized to the 
local interstellar field, Gq = Gq/I.T. Curves are shown for N22 = 1-5, Z' = 1, and two cloud 
masses M(i?co) = 3 x 10^ Mq {thick curve) and M(i?co) = 1 x 10^ Mq {thm curve). The 
radius R^,^ is where half the nuclei are in H2 (nH2/''^c = 0.25), and the radius -Rco is where 
rco = 1- For M{Rco) = 3 x 10^ Mq, Rco = 75.4 pc and for M{Rco) = 1 x 10^ M©, 
Rco — 13.8 pc. 
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Fig. 9. — Dark gas fraction /dg = [M{Rii^) — M{Rco)]/M{Rii^) versus incident radiation 
field normalized to the local interstellar field, G'q = Go/1.7. Curves are shown for N22 = 1-5, 
Z' = 1, and cloud masses M(i?co) = 3 x 10^ Mq {solid curve), 1 x 10^ Mq [dotted curve), 
3 X 10^ Mq [short-dash curve), and 1 x 10^ Mq [long-dash curve). 
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Fig. 10. — Top panel: Dark gas fraction /dg = [M{Rii^) — M{Rco)]/M{Rii2) versus incident 
radiation field normalized to the local interstellar field, G'q = Gq/1.7. Curves are shown for 
constant mean column density N22 = 1.5 and metallicities Z' = 1.9 (dashed curve), Z' — 1 
[solid curve), Z' — 0.5 [dotted curve). Bottom panel: Same as top panel for N22 — 0.75. 
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Fig. 11. — Top panel: Dark gas fraction /dg = [M{Re.^) — M{Rco)\/M{Re_^) as a function 
of the mean visual extinction through the cloud Ay = 5.26Z'N22- Curves are shown for 
constant cloud mass M{Rco) = 1 X 10^ M0, Gg = 10, and metallicities Z' = 1.9 {dashed 
curve), Z' = 0.5 {dotted curve). Clouds of higher Ay have less surface dark gas relative to 
the CO interiors. 



